Interplay of disorder and interaction in Majorana quantum wires 
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We study the interplay between disorder and interaction in one-dimensional topological super- 
conductors which carry localized Majorana zero-energy states. Using Abelian bosonization and 
the perturbative renormalization group (RG) approach, we obtain the RG-flow and the associated 
scaling dimensions of the parameters and identify the critical points of the low-energy theory. We 
predict a quantum phase transition from a topological superconducting phase to a non-topological 
localized phase, and obtain the phase boundary between these two phases as a function of the 
electron-electron interaction and the disorder strength in the nanowire. We also identify a large 
regime of stability of the topological phase in the parameter space of the model. 



PACS numbers: 71.10.Pm, 74.45+c, 74.78.Na, 74.81.-g 

Introduction. The search for topological phases of mat- 
ter has become an active and exciting pursuit in con- 
densed matter physics [Tj. Among the many impor- 
tant examples of such phases are topological supercon- 
ductors supporting zero-energy Majorana bound states 
(MBS) PH5]. A particularly promising realization of 
topological superconductivity is one-dimensional (ID) 
semiconductor/superconductor heterostructures |S]. In 
addition to being one of the simplest examples of frac- 
tionalization, zero-energy MBS quasiparticles have Ising- 
like non- Abelian braiding properties [6J and can be used 
for topological quantum computation [TJ. 

The distinct feature of topological superconductors is 
the ground-state degeneracy due to the fermion parity en- 
coded in the exponentially localized zero-energy MBS |4J. 
In a finite-length ID wire, this degeneracy is approxi- 
mate and there is an exponentially small energy splitting 
e -V£ d ue to a finite overlap of MBS |S]. Here L and £ 
are the length of the wire (spatially separating the two 
MBS localized at the edges) and superconducting coher- 
ence length, respectively. The presence of the impurities 
in the wire adversely affects the stability of the topolog- 
ical phase by introducing disorder. Indeed, each impu- 
rity in a topological p-wave superconductor leads to the 
emergence of a localized subgap Andreev bound state. 
If the concentration of such Andreev states is too large, 
the overlap between MBS changes qualitatively from ex- 
ponential to algebraic. Thus, increasing the disorder 
strength (e.g., by increasing the random impurity con- 
centration) should lead to a topological quantum phase 
transition (QPT) from the Majorana-carrying topologi- 
cal superconducting phase with quantum degeneracy to 
a trivial phase with no end-MBS in the wire |SJ[TDj. In 
addition, the question of electron-electron interaction in 
the superconducting wire may have important implica- 
tions for the topological phase, and there may be QPTs 
associated with the tuning of the interaction strength. It 
is well known that the properties of ID conductor are 



strongly affected by both electron-electron interactions 
and disorder [TT] . Clarification of their combined effect 
is crucial for our complete understanding of the topolog- 
ical phase diagram of the system and ultimately for the 
experimental realization of Majorana quantum wires in 
the laboratory, where obviously both disorder and inter- 
actions would be ever-present. 

In this Letter, we investigate an important question 
concerning the effect of both disorder and interaction on 
the stability of the topological phase and go beyond the 
non-interacting results of Refs.[9J[Tni[I2HIS]i where disor- 
der effects on the topological superconducting phase were 
considered only in the non-interacting model, and of Refs. 
|17H20| . where the effects of interaction have been studied 
in clean nanowires. We carry out a complete theoretical 
analysis in the presence of both disorder and interaction 
obtaining in the process the quantum topological phase 
diagram of the ID nanowire system. We consider a semi- 
conductor nanowire with strong spin-orbit and Zeeman 
interactions proximity-coupled to an s-wave supercon- 
ductor (SC). As shown in Ref. jS], the low-energy Hamil- 
tonian of the model in the suitable parameter regime re- 
duces to the one of an effective ID spinless p-wave SC. We 
then include the effects of both quenched disorder and in- 
teraction and, using Abelian bosonization and the replica 
method, derive a set of coupled renormalization-group 
(RG) equations for the parameters of the model. In gen- 
eral, disorder and repulsive interactions reinforce their 
detrimental effects on the topological SC phase and tend 
to eliminate the exponentially-split ground state MBS 
degeneracy associated with different fermion parity |21j . 
However, for a sufficiently strong initial induced pairing 
A, we predict a stable topological phase at low temper- 
atures, even in the presence of disorder and interaction. 
Our results shed light on the question of the stability of 
MBS in realistic situations and are important for the on- 
going experimental investigation searching for the MBS 
in semiconductor nanowire systems. 
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Theoretical model. We start with a model for a single 
channel conductor of length L with open boundary con- 
ditions, subject to Rashba SO coupling and a magnetic 
field along the wire (creating the Zeeman spin splitting) . 
It is instructive to discuss first the non-interacting theory 
for the clean system. In that case, the Hamiltonian for 
the quantum wire in the continuum is [5] 



H^ = I dx¥ 



2m 



r y d T 



V z a 2 



where h = 1, 9 = (i[>+ (x) , {x)) T is the electron spinor 
operator, a is the Rashba SO interaction parameter, <j % 
are the Pauli spin 1/2 matrices and V z is the Zeeman 
splitting. The Hamiltonian is diagonalized in the 

Fourier basis and has eigenvalues e± (k) = k 2 /2m — p, ± 



(ak) + Vg, with k the quasimomentum, and has the 
eigenmodes ^ ± (fc) defining the higher and lower disper- 
sion branches. The proximity-induced pairing interaction 
reads 



(1) 



where As is the proximity-induced s— wave gap. Our 
approach consists in expressing the pairing term in the 
v I / ± (k) basis, which leads to SC pairing interactions with 
singlet (i.e., between ^ + and subbands) and triplet 
(i.e., within the same or subband) symmetry. 
The non-trivial topological phase supporting MBS occurs 
when the condition V 2 > A| + p 2 is fulfilled [3j, which 
physically implies that only the lowest subband ^_ is 
occupied. In that case, the system is an effective realiza- 
tion of the topological phase in the form of ID spinless 
electrons with p-wave pairing 

At low energies we can linearize the band e_ (fc) around 
the Fermi points ±kp, and express the fermion field 
f_ (x) as [□] 



1 
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where the bosonic fields <f> (x) , 9 (x) are conjugate 
canonical variables obeying the commutation relation 
[4> (x) , 6 (y)] = i7rsign (y — x) /2. Physically, (f>(x) rep- 
resents slowly-varying fluctuations in the electronic den- 
sity p (x) = po — d x 4> (x) /ir, while 9 (x) is related to the 
SC order parameter. The Klein factors Ui keep track 
of the fermionic anticommutation relations and obey 

hj h U]\ = for i ^ j and UiUj = U}U l = 1. Finally, 
a ~ kp 1 is the short-distance cutoff of the continuum 
theory. The bosonized Hamiltonian Hq = Hq + Hq is 
given by 



Hq = I dx 
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Here A is the effective topological SC gap parameter, 
and for A = model Q reduces to the Luttinger liquid 
(LL) model [11] , which describes gap less plasmon excita- 
tions in the wire propagating with velocity v ~ Vp, and is 
parametrized by the the dimensionless Luttinger parame- 
ter K < 1 (K > 1) for repulsive (attractive) interactions. 
At large A, the field 9 (x) is pinned to the minima of 
sin 29 and, thus, the superconducting state breaks U(l) 
symmetry down to Z 2 . Indeed, in the limit L — > 00, 
there are two degenerate minima 9 (x) = —n/A, 37r/4, 
related to each other by the global Z2 transformation 
9 — > 9 + 7r, see Ref. [2T| for details. Such a transfor- 
mation is implemented by the fermion parity operator 



P = (-l) 



N F 



exp 



i L d x (f> (x) dx (where Np is the 



total number of fermions) , with even and odd degenerate 
eigenfunctions |e,o) = (|— 7r/4) ± |37r/4)) /y2 forming a 
fixed- fermion parity basis |21| . In Eq. ^ we have ne- 
glected the umklapp scattering which would introduce 
an additional term ~ cos (2<f) — Akpx) since we assume 
a filling incommensurate with the lattice. We are also 
taking the interaction to be short-ranged consistent with 
the canonical LL model, assuming strong screening in the 
nanowire (both by electrons in the semiconductor and 
by surrounding gate electrodes). The generalization to 
a true long-ranged unscreened Coulomb interaction does 
not change any of our conclusions as it only involves a 
very slowly-varying scale-dependent Luttinger parameter 
instead of a constant K |22| . 

We now introduce quenched disorder into model 
([2]). The Hamiltonian for short-range non-magnetic 
impurities reads = — J dx V\ (x) ^ (x) (x) (in 
the original notation), with Gaussian-distributed disor- 
der potential Vi (x) and (VJ (x) (y)) = Di,S (x — y) . 
In addition, we consider that imperfections at the 
semiconductor-superconductor interface gives rise to an 
inhomogeneous tunnel-coupling mechanism producing 
the SC proximity-effect in the nanowire [15 . In that 
case, second-order perturbation theory in the tunnel- 
coupling generates, in addition to the uniform contri- 
bution ([lj, a spatially fluctuating s-wave pairing, i.e., 
Hdis.s = ^Jdx 6 A s (x) * t (x) a v *t (x) + H.c, where 
we assume 6 Ag (x) to be a Gaussian variable obeying 
(6 As (x) 5 As (y)) — D s 8 (x — y) . Although in principle 
8 As (x) is a complex quantity encoding both amplitude 
and phase fluctuations, the latter can be neglected in the 
case of a proximate 3D SC at temperatures T <C lu p , 
where uj p — Air p^e 2 / m* is the 3D plasma frequency. 
Upon projection onto the lowest subband we obtain 
the bosonized Hamiltonians 



Hi 



dx 



-rj(x) 



V<f>{x) 
2tt 



+ £(x) 
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2na 



H.c. 



Hdis A = — / dx C (x) sin 29 (x) 
ira 



(3) 
(4) 
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Here we have defined the disordered potentials r] (x) = 
jf E ? ~o e lqX Vi (?) ,a*) = Tf E^o e lqX Vi (« " 2 M and 
C(x) = - jfY, q ~o cos {Qx)SA s (q) a/2 V z . The forward 
scattering term — 77 (x) V0 (x) /2tt can be eliminated by 
means of a gauge transformation <fr (x) — > 4> (x) — 
^ J x dy Tj(y), reflecting the fact that forward scatter- 
ing does not affect the thermodynamic properties of the 
system 125] , We next implement the replica method, 
that consists in introducing the set of "replicas" of the 
system <fr (x) , 9 (x) — ¥ <pi (x) , 6i (x), with i = 1,2, ... ,n, 
allowing a simpler integration over different disorder con- 
figurations [23] . After integrating out the Gaussian 
fields V\ and 8 A$, the replicated action becomes 




s = E 



OijSo,i — 



Db 
(2?ra) 



cos 2 [fa (x,r) - <j>j (x,r')] 



1.5 -1.0 
ln[y A (01 



0.0 



(5) 



, V {v) [ cos2[^(x,t) + ^(x,t') 

So ti = —i J drdx - (x, r) d T <j>i (x, T ) + J dT H ,i (r) . 

where the Hamiltonian Hq^ is defined in Eq. In the 
absence of SC pairing, this model has been studied in 
the context of the localization transition, predicted to 
occur at the critical value K c = 3/2, in the limit of weak 
disorder and spinless fermions [33] . For K < K c , disorder 
flows to strong coupling and the groundstate corresponds 
to a pinned charge-density- wave (PCDW) , characterized 
by a localization length £i oc oc DJ ■ Above K c , 



the LL phase remains stable, describing a "delocalized" 
electronic fluid. 

RG analysis. The critical properties of model ([5| 
can be studied in the framework of perturbative RG 
around the Luttinger liquid fixed-point. Following stan- 
dard derivations jllj . we expand the partition function 
corresponding to action S at first-order in the small pa- 
rameters Df, and D Sl and up to second order in A. We im- 
plement a RG procedure that leaves invariant the Gaus- 
sian sector of the action and obtain the following system 
of RG-flow equations 



dK\ 



dt 



dv {I 



dy A {£ 



dyb 



dl 
dy s ( I 
dl 
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Figure 1. Parametric dependence of yt {I) vs ?/a (t), as ob- 
tained from the numerical solution of the RG-flow Eqs. (JsJ) — 
(10 1, for fixed initial parameters Kq = 0.65 and y s o = 
(log-log scale). The thick dashed curve is the critical 
line, separating the topological SC phase (shaded area) from 
the non-topological disordered phase, and the thin dotted 
line is our analytical estimate yb ~ y&, valid in the limit 
{yt (£),yA (*)}-> 0. 



Here we have introduced the dimensionless variables 
y A = A/vA, y b = D b /Av 2 , y s = D s /Av 2 (with A ~ 
Uf the high-energy cutoff of the theory), and the non- 
universal numerical coefficients and B 2 of order 
unity. Physically Eq. (JgJ) describes the renormalization 
of interactions in the wire (parametrized by K {£)) in- 
duced by superconductivity and disorder. While y A (£) 
and y s {£) couple to field 6 (x) , and favor a SC ground 
state with Z 2 -symmetry, the parameter y b {£) couples to 
the dual field (f> (x) and tries to pin the density to the dis- 
order potential, thus opposing a SC ground state. These 
competing effects are reflected in the different signs of 
the prefactors in Eq. y\ (£) and y s {£) renormalize 
K {£) to larger values, inducing attractive interactions in 
the wire, while y b (£) drives K {£) — > enhancing the 
effect of repulsive interactions. We therefore see that 
repulsive interaction and disorder reinforce each other's 
detrimental effects on the topological SC. In the limit 
{yA {£) , yb {£) j Vs {(-)} — > 0, the properties of the system 
are determined by the value of K(£), i.e., the couplings 
2/A {£) and y s {£) become relevant for K (£) > 1/2 and 
K (£) > 2/3, respectively whereas y b (£) is relevant for 
K{£) < 3/2 US]. 

In order to simplify the analysis, we focus first on 
the case y s = 0. In that case, note that within the 
experimentally interesting regime 1/2 < K (£) < 3/2 
both yA (£) and y b (£) are competing perturbations flow- 
ing simultaneously to strong coupling. Moreover, in the 
non-interacting case K = 1, j/a {£) and y b {£) have the 
same scaling dimension. In order to maintain the in- 
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ternal consistency of our perturbative approach, the RG 
flow has to be stopped at a value £ max for which one 
of the couplings reaches the strong-coupling regime, i.e., 
max [y A (£ max ) , Vb (Ana*)] = 1- Although strictly speak- 
ing our approach is not applicable in the strong-coupling 
regime, the fact that 9 (x) and 4> (x) are dual fields that 
cannot order simultaneously allows us to reasonably con- 
jecture that there are no intermediate fixed-points in 
the RG flow, and therefore to classify the nature of the 
ground state according to the coupling that first reaches 
the above condition. When the two competing couplings 
reach the strong coupling regime simultaneously (i.e., 
Ub (W) = UA (4iax) = 1), the system does not order 
and this condition defines a critical line of QPTs that 
separates the topological SC phase with broken Z 2 sym- 
metry from the PCDW insulating phase (cf. dashed line 
in Fig. 0. 

From the lowest order RG equations one obtains 

Vb (£) = Vb^- 2K)l , VA (£) = yAA 2 -^)^ which to- 
gether produces the relative scaling yb ~ y^ with v = 
(3 — 2K) I (2 — K -1 ). Physically, this means that inter- 
actions (encoded in v) determine the scaling of disorder 
strength relative to the SC order parameter: for K > 1 
(attractive interactions) disorder grows slower than SC, 
while the inverse occurs for K < 1 (repulsive interac- 
tions) . In Fig. [I] we show the parametric dependence 
of yi, (£) as a function of y<\ (£), for initial parameters 
Kq = 0.65 and y s o — 0. The continuous lines corre- 
spond to the numerical solution of Eqs. (|6|-(10), and 
the dotted line is our analytical result yb ~ y A , valid in 
the limit {yb {£) , jja CO} — ► 0. At the phase boundary 
(thick dashed line in Fig. [T|), this result implies the ap- 
proximate relation ybo = y^o f° r ^ ne initial values, which 
together with the relation between impurity scattering 
time and Db'.l/r e = 2riiU^i:Afo = 2Db/vF (here m is 
the concentration of impurities, u the typical impurity 
potential and A<o = 1/ttvf the ID density of states), pro- 
duces l/2T e Ep = (A/Ep) u - Interestingly, for K = 1 we 
find that the critical condition for the topological-non- 
topological transition is l/2r e = A, which exactly co- 
incides with the results obtained in the non-interacting 

case pa EDI nun Eg. 

The results shown in Fig. [l] demonstrate that a topo- 
logical superconducting state that supports MBS could 
be in principle realized in a realistic semiconducting 
quantum wire, subject to the simultaneous effects of 
disorder and repulsive interaction, provided it is in the 
proper regime in parameter space. 

The above procedure leads to a qualitative "phase- 
diagram" in terms of the initial parameters of the model. 
In Fig. [2] we plot the critical curves in yAo-J/fco space, 
for different initial values of interaction Kq. The area 
below each curve represents the regime for which topo- 
logical SC is expected to dominate over disorder, and 
therefore stable MBS to exist. Starting from the inital 
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Figure 2. Phase diagram in yAo,Vbo space obtained for y s o = 
and different values of Ko . The curves correspond to the crit- 
ical lines yto vs j/ao, satisfying the condition yA (imax) = 
yb (fmai) = 1. The area below each curve represents the 
regime for which topological SC is expected to dominate over 
disorder. 



value Kq = 0.65, representing a strongly interacting wire, 
note that the topological region expands as the interac- 
tion becomes increasingly attractive. 

As mentioned before, for K > 3/2 and in the absence 
of pairing A = 0, disorder is irrelevant and the ground 
state of the wire corresponds to a stable LL [23] . In our 
case, we see from Eq. ^ that the LL fixed point is 
unstable against a vanishingly small pairing interaction 
yA , and directly flows to the stable topological phase [T7] . 
This result generalizes the localization transition in the 
non-SC problem to the case of topological SC. 

We now turn to the effects of a non- vanishing y s . The 
lowest order RG-flow Eq. ([6| suggests a favorable ef- 
fect on the topological SC ground state through a pos- 
itive renormalization of K (I) . Note however that this 
conclusion is only valid in the limit of small fluctuations 
Vs{£) < Va(£) (i-e., \SA s (x)\ < A s ). Actually, recent 
numerical results indicate that for typical fluctuations 
|5As(x)| larger than ~ As/4, the topological state is 
destroyed (TS]. In our description, the adverse effects on 
the topological SC are encoded in the sign-changes of the 
random potential £ (x). which induce spatial fluctuations 
in (x). From the RG flow Eq. ( 10 1 we extract a scaling 
dimension 3 — 2/ K (£) for the coupling y s (£), which for 
K (£) < 1 is smaller than that of i/a(^)- Therefore, in 



our regime of interest K (£) < 1, and for initial values 
Uso *C yAo, we expect disorder in the interface tunneling- 
amplitude to be a perturbative correction to the average 
pairing term ~ A sin (26), and our results of Eq. (|6| to be 
valid. A more general study of disorder in the interface 
tunncling-amplitude will be presented elsewhere [25 . 

Conclusions. We have carried out an RG analysis of 
the topological superconductivity in ID Majorana chain 
systems, arising in semiconductor nanowires due to the 
combined effects of spin-orbit coupling, Zccman splitting, 
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and s-wave proximate superconductivity, in the presence 
of interaction and disorder treating them on equal foot- 
ing. Our main conclusion is that the topological phase 
with non-Abelian Majorana bound states in the chain is 
stable in a large regime of the parameter space. 

We thank Chetan Nayak for valuable discussions. We 
acknowledge support from DARPA QuEST, JQI-NSF- 
PFC and Microsoft Q. 
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